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Abstract. The area of Fourier analysis connected to signal processing theory has 
undergone a rapid development in the last two decades. The aspect of this development 
that has received the most publicity is the theory of wavelets and their relatives, which 
involves expansions in terms of sets of functions generated from a single function by 
translations and dilations. However, there has also been much progress in the related 
area known as time-frequency analysis or Gabor analysis, which involves expansions 
in terms of sets of functions generated from a single function by translations and mod- 
ulations. In this area there are some questions of a concrete and practical nature whose 
study reveals connections with aspects of harmonic and functional analysis that were 
previously considered quite pure and perhaps rather exotic. In this expository paper, 
I give a survey of some of these interactions between the abstruse and the applica- 
ble. It is based on the thematic lectures which I gave at the Ninth Discussion Meet- 
ing on Harmonic Analysis at the Harish-Chandra Research Institute in Allahabad in 
October 2005. 
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1. Background and statements of theorems 

In this section we develop some basic ideas of time-frequency analysis and lead up to 
the main results to be discussed in this paper. For a more comprehensive account of this 
subject we recommend the book of Grochenig 1 14 1; see also Daubechies |5 1. 
For an L 2 function / on R d , the Fourier transform and its inversion formula 

f(co) = Jf(t) e - 27Cicot 6t, f(t) = J /(coy^'dco 

provide the expansion of / in terms of the pure 'sine waves' Q 2m(a t . (When d = 1 it is 
often appropriate to refer to t as 'time' and ft) as 'frequency.' We shall sometimes use 
these terms even when dealing with general d, for which the basic mathematical structure 
is exactly the same.) The Fourier transform is a marvelous tool, but since it involves 
the whole function / at once, it is not an efficient way to analyze the ways in which 
different frequencies enter into / at different times (as is of paramount importance, for 
example, in music). 

One way to produce a 'local' Fourier analysis of a function / is to use the windowed 
Fourier transform, which is also known as the short-time Fourier transform and is closely 
related to the cross-ambiguity function of radar theory and the Fourier-Wigner transform 
in 1101 . The idea is simple: one fixes an L 2 function on Mr such that ||0||2 = 1 and 
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considers the Fourier transform not of / but of / multiplied by translates of <j), obtaining 
a function V$f on W 1 x W 1 : 

y^/(a),x)-|/(f)^(7^)e- 2jI " 'df. (1) 

Vd,f is called the windowed Fourier transform of f with window (p. (The reasons for the 
complex conjugation on <p and the normalization || 1| 2 = 1 will become apparent shortly.) 
If one takes (j) to be a nonnegative function supported on an interval / centered at the origin 
(or at least negligibly small outside /), then Vd,f((0,x) measures how much the frequency 
co contributes to the portion of / that lives on the interval / + x centered at x. 

To obtain the inversion formula for the operator Va, we observe that the map 
(0,/) i— > V$f is the restriction to functions of the form F(x,t) — <j>(x)f(t) of the linear 
map V: L 2 (R 2d ) -> L 2 (R 2d ) defined by 

VF(co,x) = J F{t-x,t)e- 2Ki<0 'dt. 

V is the composition of the measure-preserving change of variable (x,t) i— > (f — x,t) with 
the Fourier transform in the second variable, so it is unitary on L 2 . It follows that the 
windowed Fourier transform is an isometry from L 2 (M. d ) into L 2 (R 2d ): 

||^/]] 2 =||0]] 2 ||/||2 = ]]/|| 2 . 

Therefore, we have VlV$ = I, and an easy calculation of V* then yields the inversion 
formula 

/(f) = J J V^f{co,x)${t -x)e 2KiaH dxd(0. (2) 

(As with the ordinary Fourier transform, this integral is absolutely convergent only for / 
in a dense subspace of L 2 and must be interpreted by a limiting process in general.) 

Let us look at this from another viewpoint. For co,x £ M. d we introduce the modulation 
operator M® and the translation operator T x by 

M a) f(t)=e 2 ^-'f(t), T x f(t)=f(t-x). 

Thus 

M m T x f(t)=e 2 * i '°- t f(t-x), 

T x M m f{t) = e 27Iim <'-^f(t -x) = e- 2 « ,ax M a T x f(t), (3) 
and Q and become 

Vif(m,x) = (f,M a T x 4), f = JJvif(eo,x)M a T x 4dxdm. (4) 

That is, we are using the time-frequency translates (i.e., translates and modulates) of 0, 
M a T x <j), as a 'basic set' of functions by means of which one can express an arbitrary L 2 
function /. 

However, the set {M a T x ty: co,x 6 W 1 } is highly overcomplete, and in general one can 
expand an arbitrary / using only a suitable discrete subset of it. For example, if is the 
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characteristic function of the unit cube [0, l] d , the set {MjT k (j): j,k £ 1 d } is actually an 
orthonormal basis for L 2 (R d ). (It is the basis one obtains by tiling M. d by cubes of unit 
side length with vertices in the integer lattice and using the usual Fourier basis on each 
cube.) Such discrete 'basic sets' have an obvious advantage from a computational point 
of view, as integrals must be approximated by discrete sums for numerical work anyway. 
We are thus led to the following definition. 
Given £ L 2 {M. d ) and a,/3 > 0, let 

%^,a^) = {M aj Tp k ^.j,k£ r L d }. 

3(0, a, /3) is called the Gabor system determined by (0,a,j3). (The name is in honor of 
the electrical engineer D Gabor, who suggested in his ground-breaking paper of 1 946 [13] 

2 / 2 

that S(e _7Cf ' a , 1 /cc,CC) should be a useful 'basic set'. He was not entirely correct — see 
§3.4 of 1 10 1 — but the essential idea was still a good one.) Evidently this concept can 
be generalized: one can consider {M m T x (j>: ((D,x) £ A} where A is a more general discrete 
subset of M? d . We shall consider such generalizations later, but we stick with the lattice 
A = aZ< / x/3Z< / fornow. 

The first question that must be addressed is the following: For which 0, a, and j3 does 
S (0 , a , J3 ) span L 2 (W 1 ) ? (By ' span' we mean that its finite linear span is dense in L 2 (W 1 ) .) 
This question, as it stands, is too broad to admit a reasonable answer. In particular, if one 
is given 0, the set of (a, /3) for which S(0,O!,/3) spans L 2 depends strongly on 0. For 
example, if d — 1 and (resp. 0) is supported in an interval of length I, for 3(0, a, /3) to 
span L 2 it is obviously necessary that /3 < I (resp. a < I). On the other hand, no matter 
what is, it is always necessary that the lattice al d x /3Z rf should be sufficiently dense 
in W 1 x R d . Indeed, we have the following: 

Theorem 1. If a /3 > 1, there is no £ L 2 such that 3(0, a, j3) spans L 2 (R. d ). 

When a/3 > 1 and a/3 is rational, Daubechies |4| has shown how to produce, for any 
given 0, an explicit / ^ such that / _L 3(0, OC, /3). (See also p. 107 of [5 1 for the very 
easy case where a = 1, /3 = 2.) However, when a/3 is irrational, a deeper argument is 
necessary. The reason for this rather surprising distinction lies in the structure of the group 
of operators G a p generated by the modulations and translations M a j and Tp k (j,k £ Z rf ). 
In view of (j3jl, we have 

G a .p = {^ l M a] Tp k : j,k £Z d ,l£Z}. (5) 

This group is a homomorphic image of the discrete Heisenberg group whose under- 
lying set is 1 d x Z d x Z and whose group law is 

(J, k,l)- (/, k 1 , 1') = (j + f,k+kf,l + l' + k- f) . (6) 

(This is often written with the roles of j and k switched.) Indeed, the map 

7t a ^jXl)=e- 27l,ap 'M aj Tp k (7) 

is easily seen to be a unitary representation of whose image is G a p. When a/3 is 
irrational, n a p is an isomorphism of groups. On the other hand, when a/3 is rational, say 
a/3 = p/q in lowest terms, we have 



ker(n at p)=qZ={(0,0,ql):l£Z}, 



(8) 
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and hence 

G a ,p = M d /qZ SS {(j,k,l)- j,k eZ d ,le Z/qZ), 

the group law again being given by (jfji with addition mod q in the last coordinate. Now, 
observe that Z = {(0,0,/): I G Z} is both the center and the commutator subgroup of H ( /. 
The group H</ jqZ is 'almost Abelian' : its commutator subgroup is finite, and it has normal 
Abelian subgroups of finite index (for example, {(j,qk,l): j,k G Z d , I G Z/qZ}). As a 
result, the harmonic analysis of H^/gZ can be reduced to Abelian harmonic analysis. But 
Hrf itself is a discrete group that is not 'almost Abelian,' and consequently it is not even 
type I. This means that its harmonic analysis exhibits various pathologies and involves 
queer beasts such as von Neumann factors of type II (see Chapter 7 of | llj for a fuller 
explanation of these matters). 

In §2 we develop the basic properties of the von Neumann algebra generated by G a «, 
and in §3 we discuss several proofs of Theorem 1. 

We now turn to a different question. Suppose that S(0, OC, j3) does span L 2 ; can we use 
it in an efficient way to expand an arbitrary / G L 2 1 The condition that allows everything 
to work smoothly is that S(0,a,/3) should be a frame for L 2 . The notion of frame was 
introduced by Duffin and Schaeffer in 1952, but it was not fully exploited until a third of 
a century later, when Daubechies et al 1 6 1 showed how handy frames could be in signal 
analysis. We make a brief detour into abstract Hilbert space theory to explain this idea; 
see 1 5 1 or 1 14 1 for a fuller discussion. 

A bit of notation that will be employed throughout this paper: If / is a discrete set, we 
denote functions on / by lower-case boldface letters such as c, and the value of c at i G / 
is denoted by c, . 

Let IK be a separable Hilbert space. A countable set C "K is called a frame for 

"K if there exist C\ ,C2 > such that for all / G Oi. we have the 'frame inequalities' 

Ci||/|| 2 <£|</, e ,)| 2 <C 2 ||/|| 2 . (9) 
iei 

The second inequality in (|9j means that the linear map A from "K to functions on / defined 
by {Af)i = (/, e,} is bounded from % to I 2 (I). Its adjoint A*: I 2 (I) — > Di is easily seen to be 
A*c = Y. c,ei, where the series converges unconditionally. The composition S = A* A: "K — » 
"K, given by 

is called the. frame operator. Since (Sf,f) = L|(/j<?i')| 2 > the frame inequalities (|9j are 
equivalent to the operator inequalities C\I <S< C 2 /; in particular, S is invertible, and its 
inverse satisfies C^ l I < S^ 1 < C^ l I. Since 

£K/,5- 1 e i )| 2 = £|(s- 1 /,e i -)| 2 = (5(s- 1 /),5- 1 /> = (5- 1 /,/>, 

these inequalities in turn imply that {5 _1 e,-},- e / is again a frame (with frame constants C^ 1 
and Cj~ : ); it is called the dual frame. The two frames {e,} and {5 _1 e,} can now be used 
together to produce expansion formulas for a general /G Win terms of either frame: 

/ = S(S- l f) = ^(S-'f^ei = Y,{f,S~ l ei)e u 



f = s-\sf) = s- x IMS" 1 *- 
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We now return to the Gabor system S(<j>, a, /3). If this system is a frame, an easy cal- 
culation shows that the frame operator S commutes with every M a j and 7^, so the dual 
frame is again a Gabor system, namely, S(5 _1 ^>, cc,j5); we call the dual window to 
0- 

Incidentally, Theorem 1 implies that if there exists £ L 2 (M. d ) such that S(0,a,j3) is 
a frame, then a/5 < 1, but this result is easier to prove than Theorem 1 itself (see p. 108 
of or Corollary 7.5. 1 of ltl4l l. 

In the applications of Gabor systems, one is generally interested in using windows 
with good time-frequency localization, that is, windows (j) for which both (j) and have 
reasonably rapid decay at infinity, or — what is more or less the same thing — that the 
windowed Fourier transform Vy(j)(co,x), for some nice fixed window 7, has rapid decay 
in both co and x. Experience has shown that a good quantitative measure of this decay 
is given by the norms that characterize the so-called modulation spaces M\ , which are 
defined as follows. 

A sub exponential weight on M. is a function v: M. 2d — > [0,°°) of the form v(<^) — 
q<3(p(E,)) where p is a seminorm on R 2</ and a is a nonnegative concave function on 
[0,°o) such that (7(0) =0 and \\m r ^o{r)/r = 0. (Examples: vi(«,x) = (1 + \co\ + \x\)", 
v 2 {co,x) = (1 + \co\Y', and vj,{co,x) — e^*, where a > and < b < 1.) Such a weight is 
always submultiplicative: v(^\ +£,2) < y (<=i) v (<=2)- Given a subexponential weight v, the 
modulation space M,! is defined as 

Ml = j/a 2 (l d ): |]/||i )V = jfjV y /(§)|v($)d§ <ooj, (10) 

where 7 is some fixed Schwartz-class window; M}, turns out to be independent of the 
choice of 7, as is the equivalence class of the norm ||/||i >v - Evidently the notion of mod- 
ulation space can be generalized — for example, by using an LP norm rather than the 
L 1 norm — but M\ will suffice for our purposes. For a detailed treatment of modulation 
spaces we refer to Chapter 11 of 1 14| and |8|. Incidentally, M\ (i.e., M,! where v = 1) is 
the Feichtinger algebra, often denoted by 5o(M </ ) (see |8|). 

Now, suppose is a window such that 3(0, 0£, j3) is a frame for L 2 (R d ). If (j) has good 
time-frequency localization in the sense that 6 M,i for some suitable weight v, it is 
obviously desirable that the dual window S~ 1 <p should also belong to M,l . That this is 
indeed the case is the second major result we wish to discuss. 

Theorem 2. Suppose v is a subexponential weight on M. 2d and a, j3 > 0. If € M t ! is a 
window such that S(0, OC, /3) is a frame, then the frame operator S maps M,! bijectively 
onto itself. In particular, £ M\. 

This theorem was first proved by Feichtinger and Grochenig [9 1 in the case where aj5 
is rational. The general case is a more recent result of Grochenig and Leinert 1151 : its 
proof involves some abstract machinery that was not needed for the rational case. The 
underlying reason for this dichotomy is the same as in Theorem 1 : one needs a result about 
a noncommutative convolution on l? d that is closely related to the group G a ,b defined in 
©. When aj5 is rational, these results can be obtained by Abelian harmonic analysis, but 
the general case requires a different approach. 

In more detail, the crucial ingredient for the proof of Theorem 2 is a noncommutative 
analogue of a classic result of Norbert Wiener. Wiener's theorem is usually stated as 
follows: If f is a continuous, nonvanishing, periodic function on K whose Fourier series 
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is absolutely convergent, then the Fourier series of l/f is also absolutely convergent. 
However, by passing from / to its sequence c of Fourier coefficients, this result can also 
be stated as a theorem about the convolution algebra Z 1 (Z). Indeed, taking into account the 
fact that the nonvanishing of / is equivalent to the invertibility of the operator g i— » fg on 
L 2 (R/Z), and hence of the operator a^c*aon 1 2 {T), Wiener's theorem can be restated 
as follows: Suppose c £ I 1 (Z) and the map anc*ais invertible as an operator on 1 2 {T). 
Then c is invertible in the convolution algebra l l (Z). 

The setting for the noncommutative analogue of Wiener's theorem is as follows. Given 
a real number y, we define the operation of y-twisted convolution, denoted by \\ y, on 
l\Z d x Z d ) (or l\l? d ) for short) by 

(a^b)j* = Z a lmhj-l)(k-m)^ 2Zir{J ~ l) " n - (11) 
Ijn 

We also define an involution a i— » a* r on Z 1 (I? d ) by 

( a *y) jk = a^r )e - 2 ^- k . (12) 

(Note: these definitions differ from the ones in 1 15 1 by the minus signs in the exponents. 
This is to compensate for the fact that in 1 15 1, time-frequency shifts are written as T a jM^ 
rather than M a jTp k .) I 1 (Z M ) is a Banach ^-algebra with product \ „ and involution *y. We 
shall denote this algebra by Ay: 

A y =(l X {l? A ),\\y,*y). 

Moreover, the obvious analogue of Young's inequality holds: if a S l l (Z 2d ), the operator 

L a (b) = a^ r b (13) 

is bounded on every l p (Z 2d ) with norm at most ||a||i. The analogue of Wiener's theorem 
is as follows: 

Theorem 3.1fa£ Ay, then the spectrum of the operator L a as an operator on l 2 (Z 2d ) 
is equal to its spectrum as an operator on l l (Z 2d ). In particular, if L a is invertible on 
l 2 (Z 2d ), then a is invertible in the algebra Ay. 

Section 4 is devoted to a sketch of the proofs of Theorems 2 and 3. 

2. von Neumann algebras generated by translations and modulations 

We review some notation and terminology. If "K is a Hilbert space, 23 (%) denotes the 
algebra of bounded linear operators on "K. A von Neumann algebra on "K is a weakly 
closed *-subalgebra of CB(!K). If M C 23(3i), its commutant M' is the von Neumann 
algebra of all B 6 23 (JC) that commute with all A 6 M. A fundamental theorem of von 
Neumann states that if M is itself a von Neumann algebra, then (M')' = M. Takesaki [22 1 
is a good general reference for the von Neumann algebra theory needed here. Nelson [ 18 1 
has given a particularly nice proof of the theorem just quoted. 

Given a,j3 > 0, let M a j g be the von Neumann algebra on L 2 (M. d ) generated by the 
operators and (j,k,d Z d ), that is, the von Neumann algebra generated by the 
group G a p. The first fundamental fact about M a j g is the following. 
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PROPOSITION 1. 

M 'a.p = M l/p,l/a- 

This is a special case of a theorem of Takesaki 1211 . of which there are several proofs 
in the literature (see 1201 '). It is obvious that 'M-i/p^i/a C 3Vt' a g because M m commutes 
with T x precisely when 03-xGZ.A simple proof of the reverse inclusion can be found in 
Appendix 6. 1 of Q. Given S G M.' a p and T G W-'i/p l ^ a , one shows by obtaining explicit 
representations for S and T that ST = TS; hence M' a p C (M^ 6 ^ a )' = Mi^ j/ a . 

Next, we recall that if .A is a *-subalgebra of S(JC), a faithful trace on .A is a linear 
functional T on A. such that t(AB) = t(BA) for A,Z? G .A and t(A*A) > for all nonzero 
A G A. If = 1, T is said to be normalized. 

Let ^ be the finite linear span of the operators M a fTp k ( j,k G Z d ), and define T on 

K.p^y' 

t(L c ^ M «/ 7 /30 =coo- (14) 

It is easily verified that T is a normalized faithful trace on p . Moreover, suppose 

Ri,. . . ,Rn are rectangular solids (products of intervals) in R'' whose interiors are disjoint, 
whose side lengths are all at most min(l /a,/3), and whose union is the cube [0, 1 /a] d . If 
Xn is the characteristic function of R„, then 

Af 

T(A) = a d Y,(AXn,Xn), AeM° ap . (15) 
1 

Indeed, it suffices to verify (1151 when A = MafTp^. Since the side lengths of the R„ are at 
most /3, their translates by amounts fik (k ^ 0) are disjoint, so (M a jTp k x n ,Xn) = unless 
k = 0, in which case (since |jf Rj = [0, 1 /a] d ) 

a d L{MajXn,X n ) = a d I e 2 ^'dr = S J0 . 

^ J[0.l/a] d 

Equation JT3i shows that t extends uniquely to a normalized faithful trace on M a p 
that is continuous in the weak operator topology, so that M a p is a finite von Neumann 
algebra. 

We remark that if a/3 is rational, say a/3 = p/q, then the center of M a p (that is, 
3Vt a p CiM-i/p i/ a ) is large: it contains all operators M aq jTp qk with j,k G Jf 1 . However, if 
a/3 is irrational, then the center of M a p is trivial; that is, M a p is a factor. (The ideas 
used in |7| to prove Proposition 1, as sketched above, easily yield a proof of this.) Since 
M a p has a faithful trace, it is actually a factor of type II i . 

We need one further ingredient. Given <j> G L 2 (R d ), let A^ be the map from L 2 (R d ) to 
the space of functions on l? d given by 

(A^f) jk = (f,M aj Tp^). (16) 

(We encountered such maps earlier in the discussion of frames.) If A^ is bounded from 
L 2 {R d ) to l 2 {I? d ), its adjoint is given by 

where the series is unconditionally convergent. 
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Suppose <j),Y € L 2 (M. d ) and A^,Ay are bounded from L 2 to I 2 . We can then consider 
the 'generalized frame operator' S v .<j, = A'yA ( j ) : 

S V 4f = Lif^WWocjTpkW- (17) 

An easy calculation that we leave to the reader shows that commutes with every M a j 
and Tp k , so by Proposition 1, Sy^ € 'M-i/p^/a- The explicit expansion of Sy^ in terms of 
the operators Mj/pTk/a is quite pretty; it is known as the Janssen representation: 

S v4 = (aj3)- d ^(¥,Mj/pT k/a ^)Mj/pT k /a- (18) 

j.k 

Actually, without additional hypotheses on y/ and (j) the convergence of the series on the 
right is questionable, but it is sufficient for y and (j) to belong to the Feichtinger algebra 
M\ defined by {To) with v = 1. 

PROPOSITION 2. 

If (j) and y/ are in M\, then Aa and Ay are bounded from L 2 (R d ) to l 2 (Z 2d ). Moreover, 
the Janssen representation Jl 81 for Swq is valid, and the series on the right converges 
absolutely in the operator norm. 

For the proof, see Theorem 7.2.1 and Proposition 12.1.11 of |14|. 
The algebra Mwjj i i a has a normalized faithful trace just like M. a p, which we denote 
by t'. The formula for the trace of Sy^ is very simple. 

PROPOSITION 3. 

IfAty and Ay are boundedfrom L 2 (R d ) to l 2 (Z 2d ), then 
T'{Sy 4 ) = {apT d {yA). 

This is an immediate corollary of ( fT4l i (with l/j3, 1/a in place of a, /3) and the Janssen 
representation if <f> , y g M\ . The general case is not hard to prove from formula J15l > (see 
□)• 



3. Time-frequency density of complete Gabor systems 

We recall that if M is a von Neumann algebra on J£, a vector v S "K is called cyclic for 
M if {Av: A G M} is dense in Jf, and M is called cyclic if it has a cyclic vector. The 
Gabor system S(0,a,j3) spans L 2 (M. d ) precisely when is a cyclic vector for M a | g, so 
Theorem 1 can be restated as follows: 

IfM a p is cyclic, then a[3 < 1. (19) 

The first proof of Theorem 1, in a sense, appeared before the question was even posed. 
That is, (Tf9} is a consequence of a theorem of Rieffel [ 20 1 concerning the 'coupling func- 
tion' for (a generalization of) M a p and its commutant. The connection with Rieffel's 
theorem was first pointed out in Daubechies 0. However, Rieffel's paper [ 20 1 is quite 
technical, and for those who are not specialists in operator algebras (including the present 



Some aspects of time-frequency analysis 



129 



writer) the arguments in it, and the notion of 'coupling function' itself, are hard to grasp. 
We shall not attempt to describe them further. 

At about the same time as |4|, another proof of Theorem 1, for the case d = 1, was 
given by Baggett [ 1 1. Baggett's argument is of interest because it explicitly develops the 
connection of the problem with the representation theory of the discrete Heisenberg group 
Hi. (It can be generalized to the li-dimensional case.) 

We recall that the group G a p generated by the operators M a j and Tp k is the image 
of Hi under the unitary representation % a p defined by 0. Baggett relates % a p to some 
other representations of Hi. First, let H be the Abelian subgroup of Hi consisting of 
elements of the form (j,0,l) G 1), and for 7, 8 G K let Gy g be the representation of 
Hi induced from the character % %s {j,Q,i) = e?*K8j-yl) of H Gy S acts on by 

[a^s{j 1 k,l) C ]n = ^ ni{5i - yl+n73) Cn-k- 

Next, for Tj > let Ev^ be the representation of Hi defined as the direct integral 

r® 

£ yj) = / G r , s d8, 

J[o,n) r 

which acts on L 2 ( [0, tj ) x Z) in the obvious way, and let 'Ny-q be the von Neumann algebra 
on L 2 ([0,tj) x Z) generated by the operators Ey.^/i), h G Hi. Baggett establishes the 
following facts: 

(i) n a,fi i s unitarily equivalent to E^y where 7= afi. 

(ii) Ey jT) is unitarily equivalent to Ey „/ if and only if 7 = / and 77 = rj'. 

(iii) ZNTy jj is isomorphic to 'Ny «, whenever 77 > 7. 

(iv) Ey 1 is unitarily equivalent to the representation of Hi induced from the central char- 
acter Xy(0,0,l) = e- 2ni ?<. 

The proof of Theorem 1, in the form J19I . now proceeds by contradiction. Indeed, 
suppose that 7= aj3 > 1 but M a |3 is cyclic. Then M' a ^ = Mi^ j/ a is also cyclic since 

l/a/3 < 1. (If l/a/3 < 1, it is easy to construct <j) e L 2 (R d ) such that 3(<j>, 1/j8, 1/a) is a 
frame for L 2 (M. d ); see [6 1.) Hence, by (i), 3\fy y and its commutant are both cyclic. On the 
other hand, from (iv) and the theory of induced representations, N^i and its commutant 
are also both cyclic. But then, by well-known facts about von Neumann algebras that can 
be found in |22|, the fact (iii) that Ny y and Xvi are isomorphic implies that they are 
actually unitarily equivalent, which contradicts (ii). 

A few years after Daubechies |4| and Baggett [1 1, Daubechies et al |7| found another 
proof of Theorem 1 that uses the trace t' on 'Mi/p^/a and the operators and ^ 
defined by d!6l > and dl7> in a very efficient way. (Like 1 1 1, this paper deals explicitly only 
with the case d=l, but the generalization to arbitrary d is entirely straightforward.) 

Suppose that iff G L 2 (M. d ) is a cyclic vector for M a 0. If Ay is not bounded from 
L 2 (W') tol 2 {Z 2d ), we regard it as an unbounded linear map with domain D{A Xjl ) = {/ G 
L 2 : Y.\(fiM a jTp k \i/)\ 2 < 00}; as such, it is always closed and densely defined. It fol- 
lows, by a theorem of von Neumann, that 5^.^ = A*yAy is a positive self-adjoint opera- 
tor on L 2 ^"'). As in the bounded case, it commutes with the M a j and Tjg^, and so does 
(el + Sy,y)~ l for any e > 0. It follows that if we set 



tj) = (eI + S VtV ) 
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then for any feL 2 (R d ), 

(A^f)jk = (f,M a] Tp k (el+S W)W )- l y) 

= {{eI+S w , w )- l f,Ma]T ?k w) = {A v {eI + S ¥ , ¥ )- l f) jk , 

that is, 
and hence 

S v , v (el + S v , v )- 1 =A* ¥ A$. (20) 

Next, Daubechies et al |7 1 showed that the formula in Proposition 3 for the trace of Sy^ = 
AyAfj, remains valid provided only Sy^ is a bounded operator, even if Ay or A^ is not. By 
1201 . that is the case for the y/ and <f> under consideration here, so 

(ap) d T'(S v , v (eI + S ¥ , ¥ )- 1 ) = (ap) d T'(A* v At) = ( W ,<j>). 

Now, by the spectral functional calculus, Sy.y(el + Sy.y)~ l converges strongly as e — > 
to the orthogonal projection onto the closure of the range of Sy^. But since y/ is cyclic for 
M a r, this range is dense inL 2 , so Sy^e/' + Sy.y) -1 converges strongly to I. The formula 
analogous to J15l > for %' shows that t' is continuous with respect to strong convergence, so 

(aj3)' / = (aj3)V(/) = ( V /^). 
On the other hand, we have 

<y/,<?>> = ((el+A^AjM) = e\\H 2 + \\A ¥V f > \\A v <p\\ 2 , 
and if e G l 2 (I? d ) is defined by ej k — SjoSko, 

(y/,(p) = (A* ¥ e,(p) = (e,A ¥ <p) < \\A vV \\. 

(No absolute values are needed since (y/, 0} = (ocj3 ) c/ > 0.) These two inequalities imply 
that (y/, (j)) 2 < (y, <p) and hence (aj3) d = (iff, <p) < 1. Thus the proof is complete. 

At about the same time as Daubechies et al (7), Ramanathan and Steger [19| found a 
more elementary argument to prove, and indeed generalize, Theorem 1 . They avoid von 
Neumann algebras but use an idea from a different branch of harmonic analysis: the notion 
of asymptotic density of a discrete set first exploited by Beurling in his work on balayage 
for the Fourier transform . 

To wit, let A be a discrete subset of R 2d = W 1 x W 1 , and for <E L 2 (R d ), let 30, A) be 
the corresponding Gabor system: 

S0,A)={M fl ,^: MeA}. 

We are interested in the question of whether 3(0, A) spans L 2 . For r > and G R 2d , 
let B,(r\,a) be the ball of radius r about (ri,a), and let v(r) = {2% d / d\)r 2d be its volume. 
Define 

v~(r) — min card(Anfi,-(?7,a)), 
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and define the lower density of A to be 

v-(r) 

ZT(A) =liminf— V- 

r ^°° nn 

(There is a corresponding notion of upper density, but we shall not need it.) 

We shall say that 9 ( <p , A) has the homogeneous approximation property if for every / G 
L 2 (R d ) and £ > there is an R > such that for every (rj 7 a) 6 R 2t/ there is a finite linear 
combination h of the M m T x ty's with (C0,x) G A(~)B R (ri,a) such that \\h — M n T a f\\2 < £. 
(This definition is quite a mouthful. The homogeneous approximation property implies 
that 3(0, A) spans L 2 (R d ), as one sees simply by taking (tj,a) = (0,0). But it is stronger: 
it means not only that each / G L 2 can be approximated by finite linear combinations of 
MtoTxty's with (co,x) G A, but that translates and modulates of / by arbitrary amounts 
(rj,a) can be uniformly approximated by linear combinations of M m T x ij)'s with (co,x) G A 
not too far from (77,0).) 

The main result of Ramanathan and Steger [ 19 1 is the following: 

Theorem4. IfS(<j), A) has the homogeneous approximation property, then D (A) > 1. 

Let x be the characteristic function of the unit cube. The idea of the proof is to compare 
3(0, A) with S(x,Z M ), which is an orthonormal basis for L 2 (R d ). For (77, a) G R 2d and 
r > 0, let 

K(j?,a) = linear span of {MjT k %: (j,k) G Z 2d nB r (i],a)}, 

W r {r\,a) = linear span of {Ma,T x ij): (00,x) G ADB r (r],a)}. 

Next, given R > 0, let T be the restriction to the finite-dimensional space V r (ri,a) of the 
operator Py^jj^j °Pw r+R {r},a) (Px = orthogonal projection ontoX). Since T is a composition 
of projections, its eigenvalues lie in [0,1], so its trace is dominated by its rank. Thus, 

tr(r) < dimW r+ tf(T7,a) = card(AnS r+ «(Tj,a)). 

On the other hand, the homogeneous approximation property implies that for any £ > 
we can find R > such that \\T {M jT k %) - M jT k %\\ 2 < £ for all (J,k) G I? d C\B r (r],a), 
and hence 

tr(T) = £ (TMjT kX ,MjT kX ) > (1 - e)card(Z M nB r (r],a)). 

Z 2d nB r (ri,a) 

Therefore, 

card(Z 2rf nB r (T7,a)) card(AnB r+ /j(Tj,a)) v(r + /?) 

1 v(r) " y(r + P) v^)" ' 

Pick sequences r„ — > °° and (rj n ,a n ) G K 2d such that card(AnB rn +«(TJ ra ,a n ))/v(r„ +/?)—» 
D~(A). Wehavecard(Z M nB,,,(j?„,a«))/v(r„) -> 1 and v(r„+*)/v(r„) = (r„+P) rf /r,f -> 
1 , and hence 

l-e<D"(A). 



Since £ is arbitrary, the proof is complete. 
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Ramanathan and Steger further show that if A. is a lattice and 3(0, A) spans L 2 (M. d ), 
then 3(0, A) has the homogeneous approximation property. (The proof is straightfor- 
ward.) Since ZT (A) = l/vol(M 2rf /A) in this case, we have the following generalization 
of Theorem 1 . 

Theorem 5. If A is a lattice in R 2d such that vol(R 2rf /A) > 1, there is no in L 2 (R d ) 
such that 3(0, A) spans L 2 (R d ). 

More recently, Gabardo and Han 112] have extended Theorem 1 to a very general set- 
ting using the same circle of ideas as Daubechies et al\l\, but arranged rather differently. 
They define a group-like unitary system to be a countable collection II of unitary opera- 
tors on a separable Hilbert space IK with the property that the group G(II) generated by U 
lies in {XU: U G IX, X G C}. Thus, {M a jTp k : j, k G 7r} is a group-like unitary system on 
L 2 (R d ), and any projective representation of a countable group gives rise to a group-like 
unitary system. 

Given a group-like unitary system II on JX and v G "K, one has the linear map A v from 
"K to functions on IX given by (A v w)u = (w,Uv); Gabardo and Han assume that the set 
of v G 51 such that A v is bounded from "K to l 2 (U) is dense in K. Under this condition, 
they show that 31 is the orthogonal direct sum of subspaces JX, (z G I), each of which is 
IX-invariant, and each of which possesses a vector V; such that ||w|| 2 = Y.ueU K w )^ v i)| 2 
for all w G 5X; (that is, {f/v,-: U G IX} is a 'normalized tight frame' for !H,). They show that 
Ttiei \\ v i\\ 2 depends only on IX, and they define the redundancy of IX to be 

r(U) = (Z,\\v i f)- 1 . 
They then show that: 

(i) IX has a cyclic vector (a vector v G 31 such that the linear span of {Uv: U G IX} is 
dense in 31) if and only if r(IX) > 1 . 

(ii) If the index set I is finite, then the commutant IX' is a finite von Neumann algebra. 
The formula f(A*A v ) — (w, v) (for any v and w such that A v and A w are bounded from 
JX to Z 2 (!I)) determines a faithful trace on IX', and r(IX) = f (/). 

In the case IX = {M a jTp k : j,k G Z c/ }, one can take the vectors V; to be the characteristic 
functions % n in dl5> and the subspaces 3X, to be the IX-invariant subspaces they generate. 
Then IX' = ij a , and by Proposition 3, Gabardo and Han's trace f is (afi) d %\ so 
their results (i) and (ii) imply Theorem 1 . In fact, without much additional effort, they 
yield Theorem 5 for the case of lattices of the form All 1 x B7L d where A, B G GL(n, R). 

Another generalization of Theorem 1 has been obtained by Bekka 1 2 1. Bekka considers 
an irreducible square-integrable representation n of a unimodular locally compact group 
G on a separable Hilbert space JX, and a discrete subgroup T of G such that the volume 
of G/r is finite. He proves results relating (i) the formal dimension of % (i.e., the con- 
stant d n such that ||^|| 2 ||tj|| 2 = d K f G \(E,,n(g)ri)\ 2 dg), (ii) the volume of G/r, and (iii) 
the 'center- valued von Neumann dimension' of JX, cdim(JX), as a VN(r)-module, where 
VN(r) is the von Neumann algebra on / 2 (r) generated by the left regular representation. 
(We shall not attempt to describe cdim(JX) other than to say that it is an element of the 
center of VN(r).) Bekka's abstract version of Theorem 1 in this setting is as follows: 

Theorem 6. If there is a vector v G IK such that {7t(g)v: g G T} spans JX, then 
cdim( JX) < I. 
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When specialized to the case where G is the (reduced) real Heisenberg group and T 
is the discrete Heisenberg group, this yields Theorem 5. (Note that Bekka normalizes his 
frequency variables differently, with the result that some 2k's appear in his formulas that 
are not in ours. But in fact each of his 27C's should be (2n) .) 



4. Wiener's theorem and Gabor frames 

We now turn to the proofs of Theorems 2 and 3, following the arguments of Grochenig 
and Leinert 1 15 1. We begin with Theorem 3 and then obtain Theorem 2 from it. 

Recall that we are concerned with the Banach ^-algebra Ay, which is 1 1 (Z 2d ) equipped 
with the product and involution defined by il It and (I12> . and that for a £ Ay, L a denotes 
the operator L a (b) = a \ yb. If X is a space of functions on r L 2d on which L a is a bounded 
operator, we denote the spectrum and spectral radius of L a on X by <Tj(a) and px(a), 
respectively. That is, 

°x( a ) = {A e C: A/ — L a is not invertible onX} , 

p x (a) = sup{|A|: A G ax(a)} = Jim ||L a ||^" X) . 

Here and in the sequel, 33 (X) denotes the space of bounded linear operators on X. More- 
over, we abbreviate l p (I? d ) as l p . Our goal is to prove that (7/, (a) = 0)2 (a) for all a € I 1 = 

A-y. 

Just as twisted convolution on R 2d is closely related to ordinary convolution on 
the reduced Heisenberg group (the real Heisenberg group modified so that its center is 
a circle rather than a line; see pp. 25-26 of flOll). so the twisted convolution \ y on Z 2 " 
is closely related to ordinary convolution on a modification of the discrete Heisenberg 
group H t /, namely, the group Hy whose underlying set is 7L d x 1 d x T (T = the group of 
complex numbers of modulus one) and whose group law is 

U,k,QU',k',C') = U + /,k + k>, CCe- 2 *'^'). 

(Note that (j,k,l) 1— » (j,k,e~ 2m Y ! ) is a homomorphism from H f / to Hy. It is injective if 7 
is irrational, and its kernel is qZ as in (|8jl if 7 = p/q in lowest terms.) Indeed, L 1 (Hy) is 
a Banach ^-algebra under convolution and the usual involution /*(§) = f{^~ 1 ), and the 
map J: Ay — > L 1 (Hy) defined by 

/(a)(j,fc,C) = rV (21) 

is easily seen to be a ^-isomorphism of Ay onto the subalgebra of L 1 (Hy) consisting of 
those functions / that satisfy f(j, k, Q = £~ f(j,k, 1). 

With this in mind, the ingredients from abstract harmonic analysis that are needed 
to prove Theorem 2 are as follows. In them we employ an obvious modification of the 
notation introduced earlier: if / is an element of the Banach algebra A, Oji(f) and Pa(,H 
denote its spectrum and spectral radius. The first two lemmas are theorems of Hulanicki 
1 16 1 and Ludwig IT71 . respectively. 

Lemma 1 . 1161 . Suppose § is a ^-subalgebra of a Banach ^-algebra A, and there exists a 
faithful ^-representation of A on a Hilbert space % such that if f S S and f — f* then 
=Pa(/)- Then for any f € § with f = f* we have G A (f) = C7 B(:H) (7r(/)). 
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Recall that a Banach *-algebra A is called symmetric if <7a(/) C [0,°°) for all f E A 
such that /=/*. 

Lemma 2. II17I . /f G « a locally compact nilpotent group, then L l (G) M symmetric. 

Using the fact that every locally compact nilpotent group is amenable, some standard 
facts about amenable groups, and Lemmas 1 and 2, it is not hard to deduce the next result 
(Theorem 2.8 of El). 

Lemma 3. Suppose G is a locally compact nilpotent group and f G L 1 (G) satisfies f =f*. 
Then the convolution operator Cf{g) = f * g satisfies (g))(C f) = Gn([2(G))(Cf)- 

We can now sketch the proof of Theorem 3. Suppose a G Ay. If a = a*?, we employ 
the map J defined by \2\\ to transfer the problem to L l (Hy). Since Hy is nilpotent, the 
conclusion a t i (a) = Gp. (a) follows easily from Lemma 3. The result for general a is now 
obtained by the following simple device. If L a is invertible on I 2 , then so is Lb where b = 
a*t \\y& or b = a \ y a*~i, and these b's satisfy b = b*'', so a*? \ y a and a \ Y a*v are invertible in 
I 1 by the result just proved. But then (a*? \ y a)~' \ y a*>" and a*r \\ y (a \ ya*'') -1 are left and 
right inverses for a, respectively, so a is invertible in I . The desired conclusion a)i (a) = 
0)2 (a) follows by applying this result to Ae — a where e is the identity in the algebra 
A r = l l (e jk = 8 j0 8 k0 ). 

Finally, we show how Grochenig and Leinert [ 15 1 deduce Theorem 2 from Theorem 3. 
Let us suppose that 3(0, a, j3) is a frame for L 2 (M. d ) and 6 Mj where v is a subexpo- 
nential weight function; we wish to show that the frame operator S = Sq^ (notation as in 
H7\ ) is invertible on M,! . The link with Theorem 3 comes through the following consid- 
erations. 

Lemma A. If a £ Z'(Z M ), let 

n(a)=Y,a jk M j/ pT k/o ,. 
Then % is a ^-representation of the algebra Ay where 7=1 /0!j3. 
Let 

v(j,k)=v(j/p,k/a), j,keZ d , 

and let 

ll = {ael\l 2d ):^\a jk \v(j,k)<^}. 

Since v is submultiplicative, it is easily verified that l~ is a *-subalgebra of the twisted 
convolution algebra Ay, for any y > 0. 

Lemma 5. If a £ Zp, then %{a) is bounded on M\. 

Lemma 6. If a G l\, then p t \ (a) = p t i (a). 

The proofs of Lemmas 4, 5, and 6 are all quite easy. The next one is a little deeper. 

Lemma 7. If a G l l (Z 2d ), then ll^a)!^^)) = \\L a \\ n i2 {Z 2j )y 
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The idea of the proof is as follows. Let C* (1 1 ) be the C* subalgebra of 23 (I 2 {l? d ) ) gen- 
erated by the operators L a , a G I 1 , and let C*(oc, J3) be the C* subalgebra of CB(L 2 (]R'')) 
generated by the Mj/pT^/a, j,k G Z d . One shows that the correspondence L a i— > 7f(a) 
extends to an injective *-homomorphism from C*(/ ) to C*(a,f5). But injective *- 
homomorphisms of C* algebras are always isometries. 

Lemma 8. Suppose a, j3 > 0. /fa 6 L~ ana? ^( a ) w invertible on L 2 (M ), f/ien a is invert- 
ible in the algebra l\, C Ay, w/?ere 7=1 /o:j3. 

To prove this, first suppose that a = a*?. Then 

IW*)IIb(i?) = ll L a|| S (/2) =P/2(a) = P/i(a) =Pji(a), 

where the equalities are justified by Lemma 7, the fact that a = a*?, Theorem 3, and 
Lemma 6, respectively. Thus, we can apply Hulanicki's theorem (Lemma 1) with § = 1$, 
A = Ay, % as in Lemma 4, and t K — L 2 (W l ) to conclude that a t \ (a) = cxg^j (n(a)). Since 

^ Cs(l2) (ft(a)) by assumption, a is invertible in l~. The case a 7^ a* r now follows as in 
the proof of Theorem 3. 

At last we are ready to prove Theorem 2. Recall the Janssen representation Jl 81 . which 
is valid for the frame operator S^^ under consideration by Proposition 2. In terms of the 
notation of Lemma 4, it says that 

= 7C(b), where b jk = {^,Mj/pT k/a ^>). 

Since G M,!, it is easy to verify that b G Z~, and since Q((j),a,fi) is a frame for L 2 , 
S(p,<p = ?r(b) is invertible on L 2 . Hence b is invertible in l~ by Lemma 8, so S7\ = 7r(b _1 ) 
is bounded on M\ by Lemma 5, and we are done. 

References 

[1] Baggett L, Processing a radar signal and representations of the discrete Heisenberg 

group, Colloq. Math. 60/61 (1990) 195-203 
[2] Bekka B, Square integrable representations, von Neumann algebras and an application 

to Gabor analysis, /. Fourier Anal. Appl. 10 (2004) 325-349 
[3] Beurling A, Balayage of Fourier-Stieltjes transforms, in Beurling's Collected Works, 

(Boston: Birkhauser) (1989) vol. II, pp. 341-365 
[4] Daubechies I, The wavelet transform, time-frequency localization and signal analysis, 

IEEE Trans. Inform. Theory 36 (1990) 961-1005 
[5] Daubechies I, Ten lectures on wavelets (1992) (Philadelphia: Society for Industrial and 

Applied Mathematics) 

[6] Daubechies I, Grossmann A and Meyer Y, Painless nonorthogonal expansions, Comm. 
Math. Phys. 27 (1986) 1271-1283 

[7] Daubechies I, Landau H J and Landau Z, Gabor time-frequency lattices and the Wexler- 
Raz identity, J. Fourier Anal. Appl. 1 (1995) 437^178 

[8] Feichtinger H G and Grochenig K, Gabor wavelets and the Heisenberg group: Gabor 
expansions and the short time Fourier transform from the group theoretical point of 
view, in Wavelets: A tutorial in theory and applications (ed.) C K Chui (1992) (Boston: 
Academic Press) pp. 359-397 

[9] Feichtinger H G and Grochenig K, Gabor frames and time-frequency analysis of distri- 
butions, /. Funct. Anal. 146 (1997) 464^95 



136 GBFolland 



[10] Folland G B, Harmonic analysis in phase space (1989) (Princeton, NJ: Princeton Univ. 
Press) 

[11] Folland G B, A course in abstract harmonic analysis (1995) (Boca Raton, Florida: CRC 
Press) 

[12] Gabardo J-P and Han D, Frame representations for group-like unitary operator systems, 

J. Operator Theory 49 (2003) 223-244 
[13] Gabor D, Theory of communication, J. Inst. Elec. Eng. 93(111) (1946) 429-457 
[14] Grochenig K, Foundations of time-frequency analysis (2001) (Boston: Birkhauser) 
[15] Grochenig K and Leinert M, Wiener's lemma for twisted convolution and Gabor frames, 

/. Am. Math. Soc. 17 (2004) 1-17 
[16] Hulanicki A, On the spectrum of convolution operators on groups with polynomial 

growth, Invent. Math. 17 (1972) 135-142 
[17] Ludwig J, A class of symmetric and a class of Wiener group algebras, /. Fund. Anal. 31 

(1979) 187-194 

[18] Nelson E, Notes on non-commutative integration, J. Funct. Anal. 15 (1974) 103-1 16 
[19] Ramanathan J and Steger T, Incompleteness of sparse coherent states, Appl. Comput. 

Harmon. Anal. 2 (1995) 148-153 
[20] Rieffel M, Von Neumann algebras associated with pairs of lattices in Lie groups, Math. 

Ann. 257 (1981) 403^18 
[21] Takesaki M, A generalized commutation theorem for the regular representation, Bull. 

Soc. Math. France 97 (1969) 289-297 
[22] Takesaki M, Theory of Operator Algebras I (1979) (Berlin: Springer- Verlag) 



